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Abstract: In this paper, an accurate 3-dimensional (3-D) 
analytical solution is proposed to calculate the projective 
capacitances of touch panels. In this study, both simple 
and complex patterns were investigated for the calcula- 
tion. We propose a partition strategy to divide a pattern 
into many rectangular or triangular sub-patterns. Each 
sub-pattern can be further cut into 2-D slices. The capac- 
itance of a 2-D slice was then solved by our closed-form 
formulae. The total capacitance of a pattern was obtained 
by integrating up all the partial capacitances of the slices. 
In this study, the precision of our analytical method was 
examined by comparing the simulation results obtained 
from Q3D™. The results of our method had a limited de- 
viation from the simulation results while made the com- 
putation time from one hour or some by the commercial 
software down to just a few tens of sec by Matlab™. In this 
study, the error is about 5% for simple patterns while the 
error was within 20% for a very complex pattern. 


Keywords: touch panel, capacitor pattern, analytical solu- 
tion, field-based calculation 


1 Introduction 


The touch panel (TP) has become the most popular user in- 
terface (UI) technology in modern smart products such as 
mobile phones, portable media players and personal com- 
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puters. The applications of TPs are continuously extended. 
At present, TPs have several package structures: out-cell, 
on-cell and in-cell. There are five sensing techniques such 
as resistive [1], capacitive [2, 3], optical [4], surface wave [5] 
and electromagnetic methods. In industries, the projected 
capacitive sensing is the dominant technique for commer- 
cial products. The advantages of the projected capacitive 
TP are that it can support multi-touch functions and thin- 
ner, lower-cost packages with higher sensitivity and effi- 
ciency. 

The projected capacitive TP utilizes electrical charges 
as the information carriers which are easy to interface with 
the electronic circuits. When the electrical charges on the 
sensing pads are accumulated from the driving circuit, 
the sensing circuit will discharge the electrical charges on 
the sensing pad and count how many charges by a prede- 
fined reference value. The sensitivity of the projected TP 
depends on the counts between touched and un-touched 
conditions. It is important to ensure the performance and 
signal reliability in the design phase for a good TP. Hence, 
a precise calculation for the capacitances is critically re- 
quired to predict the performance of a TP accurately before 
manufacturing. 

In industries, almost all the TP designers rely on nu- 
merical tools to calculate the capacitance of a pattern. In 
addition to conventional diamond structure, there were 
many patterns proposed to increase the sensitivity of a 
panel [6-17]. Lee et al. have analyzed the properties of a 
TP with a standard pattern of the interlocking diamond 
shape [18]. They utilized COMSOL Multiphysics™ to as- 
sess the electric fields of two conventional electrode pat- 
terns. However, numerical techniques usually consume a 
lot of computing time and memory resource even for sim- 
ple structures. Still, there are few physical meanings to 
predict possible trends or design effects from the numer- 
ical results. 

For closed-form analytical calculation for most cases, 
conformal mapping is a suitable technique for calcula- 
tions of electrostatic fields, ideal fluids, and heat flows in 
2-D problems with simple boundary conditions [19, 20]. 
Since there are few studies by closed-form analytical so- 
lution for TPs in literature, we tried to adopt some analyt- 
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ical methods which were applied for interconnects in ICs 
to this study [21-25]. 

According to the methods adopted in studies for inter- 
connections in ICs, we proposed a quasi-3D method for the 
first time by integrating the analytical solutions of the 2-D 
slices of the projected capacitive TPs. In our calculations, 
the effects of electrical field shielding and charge sharing 
were considered and integrated. 

In our previous study, our method dealt only with sym- 
metrical structures [26]. In this study, we extended the 
electrode patterns to be of arbitrary shapes. By suitable 
partition, the closed-form formulae for symmetrical struc- 
tures can be utilized for electrodes of arbitrary shapes. In 
this study, the calculations took two situations, touched 
and untouched cases, to estimate the sensitivity of a pat- 
tern. 

This paper is organized as follows: In Section 2, we will 
explain the principles of our modeling for different TP pat- 
terns. Section 3 demonstrates the results of some selected 
patterns from simple to a complex shape. In this section, 
the results will be compared with those by numerical sim- 
ulations of Q3D™. Brief discussion will be also given for 
each case. Section 4 gives the conclusion for this study. 


2 Modeling principles 


According to the development of integrated circuit (IC) 
technologies, there are a lot of 2-D analytical calculations 
for interconnects in an IC [21-27]. In our primary notion, if 
we can decompose the 3-D structure into many 2-D slices, 
we can employ these 2-D analytical methods to solve the 
distributions of the electric fields in a TP by integrating all 
the solutions of the 2-D slices. Therefore, the total capaci- 
tance of a dot on the panel can be obtained by this quasi- 
3-D method. 


2.1 Basic interconnect structures 


Let’s focus first on 2-D cases. Based on the concept of finite- 
element method, a structure of any shape can be decom- 
posed into many small pieces of simple shapes. Then, it 
can be analyzed by calculating its field equations in each 
partitioned area. By suitable partitioning, the properties of 
this structure can be approached. 

With similar concepts, the electric fields between two 
electrodes of arbitrary shapes can be estimated by decom- 
posing it into many small and simple structures. 
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Figure 1: Four basic capacitance models: (a) plate, (b) fringe, (c) 
terminal, and (d) outer plate. 


This thesis worked well in the studies by Zhao et al. 
and Zou’s group for the interconnection in ICs [27-29]. In 
their studies, they proposed that the capacitance of any 
rectangular interconnect structure can be modeled and ap- 
proximated by four types of basic linking models: plate, 
fringe, terminal, and outer plate as illustrated in Fig. 1. The 
symbols utilized in the derivation are also summarized in 
Table 1. 


Table 1: Parameters for analytical calculation. 


Symbol Parameter Definitions 

Ww ITO width 

T ITO thickness 

S ITO spacing 

Hg Glass thickness 

€0 Dielectric constant of air 
€g Dielectric constant of glass 


Fringe capacitance on the ITO 
Capacitance from ITO upper terminal 
Capacitance between parallel surfaces 
Capacitance from ITO lower terminal 
Fringe capacitance under the ITO 


Cupper-fringe 
Cupper-terminal 
Cylate 
Clower-terminal 


Ciower-fringe 


Tx The width of transmission ITO 
Rx The width of receiving ITO 
L Triangle pattern’s length 


For the plate structure in Fig. 1(a), it is composed of 
two identical conductor edges oriented in parallel, the nor- 
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malized capacitance per unit length can be obtained di- 
rectly with the width of the conductors and the spacing 
between them, i.e. 


6 T 
plate = 3 (1) 


For the fringe structure, in Fig. 1(b), of two identical con- 
ductor edges oriented vertically to each other, the nor- 
malized capacitance per unit length can be obtained by 
conformal transformation of the plate structure or sim- 
ply approximated by integrating the differential capaci- 
tance along the width of the two conductors with the cor- 
responding spacing. As proposed by Zhao et al., the inte- 
gration along simple paths would be simpler without sig- 
nificant loss in accuracy. The capacitance of this structure 
can be approximated as 


w=S+W x=S+W 


Chringe = dw _, dx 2 WwW 
E i S x ao S 


) Q 


a 
2 


w=S x=S 


The terminal structure in Fig. 1(c) counts on the relation 
between a corner and a conductor edge. The normalized 
capacitance can be derived as 


ve a x=H d 

Cterminal Ww x 4 

= ~ =In2 

é H mH+x x2 ©) 
w=0 x=0 


For the links out of the back of the smaller plate to the 
larger plane, as shown in Fig. 1(d), the capacitance would 
be 


T, 


wilh  — xlh 


dw dx 
nm) + W+x)+5 +W 
x=0 


Couter-plate - 


4 nT (4) 
2(1 + mM +W) 


Based on the line integrating results of the above four 
basic structures, the capacitances between conductors in 
rectangular structures can be estimated by the related for- 
mulae. 

By the modeling proposed in refs. [27-29], the ca- 
pacitances of several regular interconnect structures in 
nanometer-scaled ICs can be approximated by summing 
up the partial capacitance of each basic structure. Accord- 
ing to the similarity between the metal interconnection in 
ICs and in TPs, it is believed that the partition method for 
the interconnect in ICs can also be applied on those of TPs. 
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2.2 Cutting electrode structures into 2-D 
slices 


In today’s TP technology, a general TP structure is com- 
posed of two parts: the indium tin oxide (ITO) traces in 
X and Y directions. The cross-over of the two directional 
traces makes an array of capacitors to sense the variations 
of charges around a local area on the panel. Let’s take the 
popular diamond structure in Fig. 1 which was proposed 
by Apple Inc. in 2013 as the first example. As seen in Fig. 2, 
the red diamond electrodes are connected horizontally in 
X ITO traces and the blue ones are in vertical Y traces. For 
most commercial products, the ITO traces in both direc- 
tions are grown at the same layer to form a coplanar struc- 
ture. In order to prevent possible shorts between X and Y 
traces, bridges are introduced as shown in Fig. 3(a). Since 
each crossover point has four neighboring diamonds and 
this structure is repeated along the panel, these diamonds, 
as marked by a circle in Fig. 2, can be a basic unit in calcu- 
lation with the repeated boundary conditions. 


Figure 2: Top view of a diamond pattern for capacitive TPs. 


Figure 3(b) depicts the cross-sectional structure along 
the cutting line in Fig. 3(a). We can make many slices by 
cutting lines on the surface in Fig. 3(a). Figure 4 illustrates 
our method to produce many slices of symmetrical struc- 
tures. Since the electrode pattern is symmetrical in Fig. 4, 
we can divide this structure into the same 4 parts of the 
same capacitance. This step will make the calculation sim- 
pler. In the meanwhile, it should be noted that the effect of 
bridge was neglected in this study. For the triangular ABC 
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Figure 3: (a) Perspective view of the structure of the sensing pat- 
tern at the center of the marked circle in Fig. 2. A bridge is used to 
connect the pads of the same trace. (b) Cross-sectional view for the 
layer structure along the cutting line in (a). The insulate layer under 
the bridge is usually the same as the cover glass. 


in Fig. 4, we can cut further to make the X- and Y-ITO in 
each slice have the same geometry. However, it should be 
noted that the width of the conductors in every slice, Wi, 
will change from one slice to another. 

For each cutting line, the related cross-sectional struc- 
ture will look like that in Fig. 5. After the capacitance of 
a 2-D slice in Fig. 5 is obtained, the total capacitance in 
Fig. 4 can be obtained by integrating all partial slice capac- 
itances. The 2-D slice as shown in Fig. 5 was the structure 
basis in our calculations. We tried to partition all electrode 
patterns in this study and cut them into the basic slice as 
in Fig. 5 with different geometry parameters. 


2.3 Calculation for 2-D symmetrical 
structures 


For the basic 2-D slice as shown in Fig. 5 and Fig. 6(a), there 
are 6 parts to be calculated. Since it is still a symmetric 
structure, it can be folded by regarding there is a vertical 
conducting plan in the center. We can calculate the partial 
capacitances for a half part as shown in Fig. 6(b). There- 
fore, C, in Fig. 5 is 1/2 of C upper-fringe in Fig. 6(b). In this 
figure, one can recognize C is equal to a serial connection 
of two fringe structures as in Fig. 1(b). Therefore, we can 
denote the related capacitance as Cypper-fringe aS Show in 
Fig. 6(b). The other partial capacitances in Fig. 6(b) can be 
analyzed in the same way. All the capacitance values from 
C, to Ce can be derived from the corresponding results in 
Fig. 6(b) by half. 

It should be noted that there might be different di- 
electrics in an area which will make the calculation a little 
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Figure 4: Cutting for symmetrical slices in triangle ABC for the dia- 
mond pattern. 
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Figure 5: Extended cross-sectional view of the basic slice marked by 
the dashed rectangle in Fig. 3. This structure can be divided into 6 
parts. The geometry parameters are also denoted. 


complex. For example, the area for C,in Fig. 5 may contain 
cover glass and air if the thickness of the cover glass, Hg, 
is less than the width of the electrode, W. Therefore, eqs. 
(1)-(4) are too simple to be applied for TP structures. The 
problems will be solved in the next sections. 


2.4 Partial capacitances for the basic slice 


For the structure in Fig. 5 or Fig. 6(b), let’s assume all the 
coplanar conductors are surrounded by the same dielec- 
tric material, for example €gjejectric- According to the parti- 
tion in Fig. 6(b), the total capacitance C for the basic struc- 


Unauthenticated 
Download Date | 7/20/17 3:06 PM 


DE GRUYTER OPEN 


——— 


en” oi oe oe 


: ‘“ “NX . 
CTyinge —* \ Cupper-fringe * \o ‘outer-termina\ 


/ = \ = \ ! 
Lf CThermind Cup)er-terminat | A 
. . . _ 


Ghat rT Cc 
s/2 outer-plate 


CT puter-plate 
ie 
Croder-terminat —! 


\ 
—_— kp 7 \ 


Ciower-fringe , t Cc 


‘sitet teraitiall 


Figure 6: (a) Dimensions and partial capacitances for the basic slice 
in Fig. 5 with two identical conductors in parallel. (b) Partitions for 
the half structure by Zhao’s models. The adopted models are also 
indicated in the subscripts. 


ture in Fig. 5 can be derived by 
C a 4 
2 E dielectric 


E dielectric 


Cupper-fringe + Cupper-terminal + Cylate + Clower-terminal 
+Ciower-fringe + Couter-plate a 2 C outer-terminal 


(5) 
It should be noticed that Cypper-terminal cannot be ob- 
tained simply by eq. (3) due to its different linked area. 
To account for the charge distribution and to compensate 
the field distortion due to adjacent plate and fringe capac- 
itances, Zhao et al. [27, 28] proposed an empirical approx- 
imation for these two corners, i.e. 


Ciower-ter minal ~ 2 
E 


(6) 


and 
Cupper-ter minal ~~ 1 


E 70 (7) 


Therefore, the total capacitance for the basic slice in Fig. 5 
or Fig. 6(a) is 


C T 2 ( x) 3 
—— =5+—=In/{1+—]+= 
S a7 


€ dielectric 


1 nT 
a Gerecent (8) 


For the TP structures, however, the dielectric layers 
are different above and below the ITO conductors. The 
above formulae should be modified. In our calculation, 
the dielectric constant was taken as a variable along the 
integrating path in eqs.(1)-(4). The derived partial capac- 
itances for the basic slice in Fig. 6(b) are summarized in 
Table 2 [26, 30]. For present study, the interface effect on 
the flow lines of electric fields was neglected. In Table 2, 
Eoft-fringe ANG E er termian ate effective dielectric constants 
by considering the field flow lines passing through glass 
and the air. 
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2.5 Calculation for symmetrical structures 


We calculated for two situations: untouched and touched. 
In general, the dimensions of a capacitive sensing pattern 
need to take that of a human’s finger into consideration. 
For a too large capacitive sensor on the panel, the user’s 
finger will touch over less than a dot unit. In this case, the 
sensitivity and linearity of the TP would be deteriorated. 
On the other hand, it is not cost effective to make the dots 
on the panel too small. For most commercial products, the 
dimensions of a capacitive sensor dot are designed to be at 
least 50% of the width of a finger. Therefore, in this study, 
we assumed the basic slice was in a size that can be cov- 
ered by one finger only. In our calculations, a finger on the 
panel was represented by a grounded post to be placed on 
the cover glass. 

The basic structure in our calculation can be repre- 
sented in Fig. 7. Fig. 7(a) depicts the schematic structure 
for the untouched situation. Fig. 7(b) emulates the touched 
situation. The finger was modeled as a grounded metallic 
post on the surface. Due to the redistribution of electric 
fields under the finger, the integration range will be con- 
fined under the finger. Therefore, the partitions for the ba- 
sic structure under untouched and touched situations will 
be different. 


Figure 7: Schematic structures for (a) untouched and (b) touched 
situations. 


Figure 8 illustrates more detailed partition and shows 
the scenario of the touched case. In this case, some cou- 
pling effects become prominent. When the finger touches 
on the cover glass, the distribution of the electrical fields 
will change. Figure 9 depicts the effects of field shielding 
and charge sharing as the following. 

The first is field shielding as illustrated in Fig. 9(a) 
which the effective coupling area between the two ITO sur- 
faces is reduced by the approaching of an additional finger. 
The second is charge sharing as in Fig. 9(b). When there 
are multiple conductors, the flow lines of the electric fields 
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Table 2: Closed-form formulae for partial capacitances in Fig. 6(b). 


Partial 
Capacitance 


For uniform e€ 


and glass €g 


For untouched TPs with air € 


For touched TPs with air €9 and 
glass €g anda pod 


€o€gzIn(1+ 2W/S), for Hg < S/2+W 


H,-S/2 
S72 


€o€gt In (1 + 
Cupper-fringe : In (1 + wn) 1 : 
+7 eff-fringeEO In(1 + 


for Hg > S/2+W 


W+Hy+S/2) 
S/2 


1 (S + 2W1)(S + 2Hg/pi) 
€0€g (4 In ( S(S + 2H¢/pi + 2W;) )) 


we €0€eq, for Hg 2 S/2 (2 In(1 + 1.297H,¢1/S))? 
Cupper-terminal aed Eo€g 5 rm 
€ef-terminal€o x» for Hg < S/2 <In(1 + 1.2974H,)/S) + 4 1n(1 + 0.32445/Hg) 
T T T_2 Hg+25/n+T- 

Ci ; ~ 2€ £0 €0 

lower-terminal ra Tt id 

Clower-fringe £In(1 + 372) £in(1 + 372) £1n(1 + 372) 

a €ot €eff-termina T T 

Couter-plate 2 ea 4 (Sema) In (1+ sesahtgreamy) 941m (1+ zasaftereem) 
Couter-terminal = & = £ 2 € 

Note: 


1. €og-fringe = 1+ wey 2. €og-terminal = 1+ ane 3. Hg: = min(Hg,S/2) 4.1, = min(T, ,/S? + H2 - Hg) 


min(W, Hg-S/2), for Hg =S/2 
W,= 
O for Hg < S/2 


Figure 8: Partition for partial capacitances for the touched situation. 


wll 


(a) (b) 


Figure 9: The integration paths change due to the effect of (a) field 
shielding and (b) charge sharing. 


Vv Ss x 


Figure 10: Another basic structure with asymmetrical trace widths of 
X and. 


may diverge to neighboring conductor surfaces rather than 
a single surface. 

In Fig. 8, there are five parts to be calculated. For con- 
venient calculation, the finger was presumed as an infinite 
ground plane to simplify the mathematical forms of the 
partial capacitances. Therefore, the integration for every 
partial capacitance needed to be recalculated to count the 


effect of the finger. The derived formulae are summarized 
in Table 2. 
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(b) 


Figure 11: Slicing strategies for (a) diamond pattern, (b) triangle 
pattern, and (c) concentric square pattern. 


2.6 Calculation for asymmetrical structures 


The above discussion deals only for the basic structure of 
two ITO traces of the same widths. By suitable transfor- 
mation, our formulae in Table 2 can also be applied for 
asymmetrical structures. For the asymmetrical structure 
as shown in Fig. 10 with unequal widths, we can employ 
suitable transformation of the unequal-width structure to 
utilize the formulae in Table 2 by an effective width [27]. 
And the effective width, W, in Table 2 should be expressed 


as 
(\ /XY(X + S\(¥ +S) + xY) 
(X+Y+S) ; 
where X and Y represent the widths of the two ITO pads in 
a slice. 

Therefore, with the transformation by eq. (9), we can 
also utilize the same formulae in Table 2 by making slices 
of asymmetrical structures for electrode patterns. The ba- 
sic structure and the formulae in Table 2 were used for cal- 
culation of the patterns in our study. 


(9) 


2.7 Slicing/Integration strategies 


In this study, five patterns from simple to complex were 
chosen for the investigation. Figure 11 depicts the slicing 
strategies for 3 simple patterns to make sliced pieces have 
a basic structure as shown in Fig. 5 or Fig. 10. For each slice 
with its capacitance per unit length, the total capacitance 
for an electrode pattern in Fig. 11 can be integrated as the 
following. 

For the diamond pattern in Fig. 11(a), the slice capaci- 
tance can be denoted as C(W)), its total capacitance can be 
calculated as 


N 
Croat = [ COW)AL™ Y>CCW)AL —— (00) 
it 


For the triangle pattern as show in Fig. 11(b), the total ca- 
pacitance can be approximated as 


N 
Croat = f COX VAL YOK, VAL — (A 


i=1 


With the transformation by eq. (9), one can use the formu- 
lae in Table 2 easily to obtain the total capacitance. 

The slicing strategy was a little complex for the con- 
centric square pattern proposed by Yu et al. [9] in Fig. 11(c). 
Figure 12(a) shows our strategy by dividing the pattern 
into two parts. For part I, we cut this part into asymmet- 
ric slices, as depicted in Fig. 12(b). 
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For part II, the sliced pieces are triangles denoted as 
the original slice in Fig. 12(c). To calculate its capacitance, 
the direct way is to recalculate the line integrations for the 
related partial capacitances in a way similar to those in 
eqs.(1)-(4) [31]. We can transform the original slice to a tri- 
angular slice as the reformed slice in Fig. 12(c) by presum- 
ing the capacitances of both slice are very close. In addi- 
tion, if the lengths X; and Y; in Fig. 12(c) are similar, the 
capacitance of the reformed slice can be approximated by 
half of the capacitance of the basic slice as the slice III in 
Fig. 12(c). Therefore, by suitable partition without losing 
too much accuracy, the original triangular slice in Fig. 12(c) 
can be reformed by regarding it as a symmetrical triangu- 
lar slice. Therefore, the capacitance of slice II C;_;; is pre- 
sumed to be close to C;. And C; 4, is roughly half of the ca- 
pacitance of the basic structure, slice III. 

In our study, the number of triangular slices was not 
many. The accumulated error in the calculation for trian- 
gular parts can be controlled in a limit. 

With the basic slice, symmetric or asymmetric, and the 
triangular slice, we can divide an electrode pattern into a 
combination of rectangular and triangular patterns. The 
partition strategy for complex patterns will be shown di- 
rectly in the next section. 


3 Results and Discussion 


In this study, five TP electrode patterns were selected to 
investigate the precision of our analytical calculation. In 
addition to the mentioned diamond, triangle, and concen- 
tric square patterns in previous section, two special pat- 
terns proposed by HTC Inc. [10] and Raydium Semiconduc- 
tor Inc. [12, 14] were also investigated. 

Our calculation was performed by Matlab™. The pre- 
cision of our method was verified by comparing the simu- 
lation results by Q3D™ which is popular in industries for 
estimating the capacitance of TPs. 

In this study, the capacitance of the bridge in each pat- 
tern was not considered due to its little effect on the sen- 
sitivity, AC. The sensitivity is defined as the difference be- 
tween the capacitance without touch, Cy7, and that under 
touch, Cr. 


3.1 Diamond pattern 


As indicated in Fig. 11(a), the slicing strategy makes each 
slice having the symmetrical structure as in Fig. 6. Fig- 
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F f 
f IT. 


I: Original IL: Reformed Ill: Basic 
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Figure 12: Slicing strategies for concentric square pattern. (a) The 
pattern is divided into two parts. (b) Slicing for part I. (c) Slicing 
for part Il. The original slice (I) can be approximated by reforming 
it to a symmetrical triangular shape (II). This reformed slice can be 
estimated by calculating the capacitance of slice Ill. If X ~ Yj, the 
formulae in Table 2 can be applied. 


ure 13 also shows five structural variables that can be ad- 
justed, ie. W, S, T, Hg and eg. 

In this calculation, S was less than 100 nm for optical 
requirements. T was less than 1 pm for ITO layers. For all 
the cases in this study, Hg was selected as 0.7 mm. And €g 
was selected as 7 for tempered cover glass. 

Fig. 14 demonstrates the comparison of the analytical 
calculation with the numerical simulation by Q3D™. As 
seen in this figure, both the analytical and numerical re- 
sults are very close. In this calculation, the maximal de- 
viation between two results was within 3%. Comparison 


Unauthenticated 
Download Date | 7/20/17 3:06 PM 


DE GRUYTER OPEN 


=e 
_--" ~~ 


Glass ¢ = eg 
ITO and Finger= perfect conductor 


Figure 13: Structural parameters for the diamond pattern. 
, Capacitance value of the Variable | 


S —=— Simulation by Q3p™ 
—© = Calculation by Matlab™ 


Figure 14: Comparison of Cyr, Cr, and AC for the diamond pattern 
as a function of spacing, S. The analytical calculation results are 
also compared with those simulation results by Q3D™. In this case, 
W=5mm,T=1pm. 


for other parameters was also studied in our previous re- 
port [26]. 


3.2 Triangular pattern 


Figure 11(b) indicates the strategy of slicing for the triangu- 
lar pattern. Figure 15 shows the structural parameters that 
can be changed in this investigation. 

The selected parameters are similar to those for the 
diamond pattern and indicated in Fig. 16. While estimat- 
ing the effects of other parameters on the capacitances, 
W =5 mm, and L = 20 mm. As seen in Fig. 16, the devia- 
tion between our analytical calculation and the numerical 
results is within 5%. In this case, L was fixed as 20 mm 
while W was changed from 2.5 mm down to 1 mm. It can 
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Figure 15: Structural parameters for the triangular pattern. 


Figure 16: Effect of dimension ratio L/W on Cyr, Cr, and AC for 
the triangular pattern. In this case, L = 20 mm, S = 30 um, and 
T=1pm. 


be confirmed that the transformation by eq. (9) is good for 
the calculation for the triangular patterns. 


3.3 Concentric square pattern[9] 


The structural parameters for the concentric square pat- 
tern were indicated in Fig. 17 for this calculation. Figure 18 
is the effect of dimensional ratio of Tx/Rx where we kept 
the total length W =5 mm. In this figure, it should be noted 
that the maximal deviation between the analytical results 
and the simulation results is 5%. 

As seen in Fig. 18, the worst case may be found for 
small Tx/Rx ratio. The precision of the analytical calcu- 
lation is good for Tx/Rx larger than 4. In this figure, the 
capacitances of Cy and Cyr are underestimated as com- 
pared to those by numerical calculation. Even our calcu- 
lation may arise more errors, it seems that the sensitivity 
AC can tolerate smaller Tx/Rx ratio down to 1 or 2. This 
finding indicates our approximation in Fig. 12 would re- 
sult in larger errors in the calculation. One should be care- 
ful while making a triangular sub-pattern as in Fig. 12. For 
partitions keeping Tx/Rx higher than 2, the precision of the 
analytical calculation can still be acceptable. 
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Figure 17: Structural parameters for the concentric square pattern. 
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Figure 18: Effect of dimension ratio Tx/Rx on Cyr, Cr, and AC for 
the concentric square pattern. In this case, W=5 mm, S=100 um, 
T=0.1 pm. 


3.4 Hexagonal pattern[10] 


This pattern was proposed by HTC Inc. in 2013. In Fig. 19, 
the structural parameters and the partition strategy are in- 
dicated. As indicated in this figure, by suitable partition, 
this pattern can be divided into many rectangular and tri- 
angular sub-patterns. By this partition strategy, the formu- 
lae in Table 2 can be employed for these sub-patterns to 
calculate the total capacitance. Figure 20 shows the effect 
of dimensional ratio of Tx/Rx. In this figure, it should be 
noted that the maximal deviation between the analytical 
results and the simulation results is 5%. 

As seen in Fig. 20, similar results can be found as in 
Fig. 18. Capacitances Cy and Cyr by our calculation are 
a little lower than those by numerical calculation. How- 
ever, the sensitivity AC can tolerate a wider range of geom- 
etry ratios. While checking the partition in Fig. 19, we may 
attribute that most of the errors come from the triangular 
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sub-patterns. In this case, the approximation in Fig. 12 still 
works fine with a 5% deviation from the numerical calcu- 
lation. 


3.5 Complex pattern[Y07,Y09] 


In addition to HTC’s pattern, we also investigate a rather 
complex pattern proposed by Raydium Semiconductor 
Inc. In Fig. 21, the structural parameters and the partition 
strategy are indicated. Since this pattern is complex, the 
partition lines in this figure divide this pattern into many 
fragments or sub-patterns. Figure 22 compares the results 
from the analytical calculation and the numerical simula- 
tion as a function of the length of the hand. In this figure, 
the largest deviation between the results by the two tools 
was around within 20%. 

The larger error may be attributed from the calculation 
of triangular sub-patterns. As checking our partition in 
Fig. 21, one can find there are many triangular sub-patterns 
with a smaller geometrical ratio for Tx/Rx in Fig. 18. These 
triangular sub-patterns may deteriorate the resulted preci- 
sion of the calculation. 

This finding may indicate that our partition strategy is 
not very efficient for this complex pattern even our analyt- 
ical results are still close to those by numerical simulation. 
If we can reduce the number of triangular sub-patterns, it 
is believed the precision of our analytical calculation can 
be improved. Nevertheless, our findings from the results 
of different patterns would convince us that our analytical 
calculation has a precision competing to commercial nu- 
merical tools. 


Tx+Rx=900um. 


Figure 19: Structural parameters for the hexagonal pattern. The 
dashed line indicates how this pattern was divided and calculated 
with different slice structures. 
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Figure 20: Effect of dimension ratio Tx/Rx on Cyr, Cr, and AC for 
the hexagonal pattern. In this case, Tx+Rx= 0.9 mm. The other pa- 
rameters can be found in Fig. 19. 


3.6 Discussion 


According to the above results in Figs. 12-22, it can be 
found that our calculations were accurate enough and can 
compete to commercial software tools which utilizing nu- 
merical techniques. 

With our analytical formulation and partition strate- 
gies, it is convenient for designers to evaluate the quan- 
titative properties and the effects of design parameters of 
a TP pattern. Especially, the computing time was largely 
reduced from several hours to some tens of seconds for 
the same structure. Since the analytical method costs lit- 
tle computational resource, this method is suitable for 
quick evaluation of effects of design parameters in the 
early stages of designing an ITO pattern. In other words, 
an optimal set of factors can be obtained by our analytical 
method. 


4 Conclusion 


This paper presents a field-based model with partition 
strategies to calculate the capacitances of TPs from sim- 
ple patterns to a complex structure. The formulae for ba- 
sic slices were derived. For complex patterns, one should 
be careful to make triangular sub-patterns to control the 
precision of the analytical calculation. By comparing with 
the results of Q3D™ simulations for the same TP pattern, 
our analytical model demonstrated an excellent accuracy 
for the capacitance calculation. 
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Figure 21: Structural parameters for the selected complex pattern. 
The dashed line indicates our partition strategy to divide this pat- 
tern into many triangular and rectangular sub-patterns. 
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Figure 22: Effect of the length of the hand in Fig. 20 on Cy7, Cr, 
and AC for Raydium complex pattern. In this case, the employed 
parameters are noted in Fig. 21. 
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